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Proofs are developed to explicitly show that the ionization energy theory (lET) is a renormalized 

theory, which mathematically exactly satisfies the renormalization group formalisms developed 

by Gell-Mann-Low, Shankar and Zinn- Justin. However, the cutoff parameter for lET relies on 

the energy-level spacing, instead of lattice spacing in k-space. Subsequently, we apply the earlier 

proofs to prove that the mathematical structure of the ionization-energy dressed electron-electron 

screened Coulomb potential is exactly the same as the ionization-energy dressed electron-phonon 

("V^ ' interaction potential. The latter proof is proven by means of the second-order time-independent 

Cn \ perturbation theory with the heavier effective mass condition, as required by the electron-electron 

►^.,^, screened Coulomb potential. The outcome of this proof is that we can derive the heat capacity and 

j^ ' the Debye frequency as a function of ionization energy, which can be applied in strongly correlated 

matter and nanostructures. 
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i. 

C/5 , I. INTRODUCTION 

The specific application of the renormalization group theory (RGT) in high energy physics requires removing the 
ultraviolet (fc — ^ oo) divergences, where fc is a wavenumbcr. Whereas, in condensed matter physics, RGT applies to 
i-rt . the study of phase transitions of spin-dependent (magnetic) system and removing infrared (fc — >■ 0) divergences. For 
r^ ' example [l| , in a d-dimensional integral, 

O' 
Q 

r^ dfik _ f^ k'^-^dk J 

^ ■ X k^ + m^ ~ Jo k^+m^ "" 1 ,a 2 ' ^' 

>- ■ Ji) < J[) < l^jj^^ j^^ ^^2 and 4;i^ for d<2 

in 

''^ ' where A is a cutoff parameter in a reciprocal space (tt/o). For large K/m? and d > 2, Eq. ([T]) diverges like A''"^ 

in the limit, A — ?> oo (ultraviolet divergence). Here, m denotes mass and note that it is common in the high-energy 
physics literature to take h ~ c = 1, where H and c are the planck constant divided by 27r and the speed of light in 

C~^ ' vacuo, respectively. Therefore, fc^ -I- m^ = c^fi^k^ + m^c^ = E^ . On the other hand, for d < 2 and to = 0, the integral 
diverges like In A (d = 2) and ^_2 (d < 2) as A — > (infrared divergence). Hence, it is appropriate to start with a 

f~*) \ general formalism of RGT (given below) based on the work of Zinn- Justin Q and Goldenfeld [3[ . From this general 
approach, we will be able to identify the relation between lET and RGT, which will be discussed in the subsequent 
sections. 

><■ 

5t . A. General formulation: Zinn-Justin formalism 

The renormalization procedure originated as a an empirical program to overcome the problem of infinities (singu- 
larities) in quantum electrodynamics [2| . The renormalization procedures have enabled one to calculate finite values 
from divergent expressions. In addition, the predicted finite values obtained from these renormalized results have 
been confirmed by experiments with high precision, which eventually gave birth to the formal RGT as it is known 
today y. In the subsequent sections, the RGT is introduced by means of the Hamiltonian flow. This means that 
the original Hamiltonian is renormalized with an arbitrary cutoff parameter, A, and then, in the flnal formula, A is 
removed via the limits, A — >■ or A — > oo. All these will be clarified shortly. 
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1. Hamiltonian flow 

We first introduce tlie parameter, A, wliich is the dilatation parameter in such a way that the effective Hamiltonian 
is given by H{^. Therefore, H i— >■ H/^ where i^A is the renormalization group (RG) transformed Hamiltonian Q] and 
H^ denotes mapping or creating (within the Hamiltonian space) a function Hf^ from H that depends on A explicitly. 
In this case, the RG idea is that the bare Hamiltonian, H = H\=i = Hi, while the RG transformed Hamiltonian 
(i?A) and for A > 0, we need to know the flow of the H/^. The general form of the flow equation is given by [2] 

A-^HA = Tm], (2) 

where 7" is the transformation in the space of the Hamiltonian. We assume [4I that the transformation is Markovian 
in which, the transformation does not depend on the trajectory between H and H\. In addition, the Markovian process 
is stationary, meaning that the right-hand side (RHS) of Eq. ([2]) is not an explicit function of A, but only through 
H\. Moreover, the identity, 

aA ' 



dA dluA' 

will be used without warning from now onwards. Here, Eq. ([5]) is a RG equation in which A determines the scale 
changes in the dynamical process. A fixed point can be defined as independent of A where we label H* as the 
Hamiltonian at a fixed point, in which Eq. (j2]) can be written as [21 

A-^Hl = T[HX] = 0. (3) 

Let us now work near a fixed point, in other words, we can perform linearization: H\ = H* + AH\. Linearizing 
Eq. ^ leads us to [H 

A-^AHa = C*{AHa), (4) 

where C* is a linear operator independent of A. The solution to Eq. Q is Q] 

Ai/A = A^*Ai/A=i, (5) 

from Eq. ([S]) , it is possible to study the RG flow near the fixed point where Eq. ([5]) gives us the relationship between 
the transformed Hamiltonian, AH/^^i somewhere in the neighborhood of H* and AHi, also in the neighborhood of 
H*. The fixed points are also known as the critical (fixed) points where the fixed-point Hamiltonian (H*) is invariant 
to A. For example, there are two trivial critical points for a ferromagnet, for T = (all spins are aligned) and for 
T — >■ 00 (all spins are randomly oriented). The non-trivial critical point is for T = Tc where, below Tc, some spins 
are aligned that gives rise to the correlation length. The correlation length is defined as the length for a group of 
spins that are aligned below Tc- Any RG transformations must give same results at T = (ordered) and for T ^> c« 
(disordered). 

Note 1: Every solution obtained from the lET Hamiltonian, is a renormalized solution (can be proven using the 
Zinn- Justin general formalism given above.) 

2. Classification of eigenvectors 

Assuming [2| that C* has discreet spectra with eigenvalues, h, which corresponds to a set of eigenvectors, 3^^. 
Therefore, one can expand Ai^A on the eigenvectors, 3^^ of £*: 

Ai7A = ^/i»(A)X. (6) 



The flow equation, Eq. ^ is now given by [2| 



A^/i.(A) = khdA), (7) 



which has the solution, 



h,{A)^A''h,il), (8) 

Here, the eigenvectors, 3^,; can be classified according to the sign on the eigenvalues, k, provided that A > 0, which 
is true in our case because, 1 < A < oo. Note here that we deliberately chose the lower limit A = 1, and not A = as 
mentioned previously. This choice does not affect the physical theory in any way since we can set any arbitrary values 
for A, which are convenient to us so as to capture the required physics accurately. Furthermore, in this case, setting 
A = is meaningless as can be seen from Eq. (|S]). Having explained that, now we can return to the classification 
issue, which are given by: 

(i) li > 0: hi and the components of hi will grow or increase with A, and the system moves away from the fixed 
point Hamiltonian, H* (repelling or unstable fixed point). In this case, hi is said to be relevant because increasing 
A (increasing dilatation) decreases effective correlation length. For example, the Curie temperature of a ferromagnet 
is an unstable critical fixed point, which means that any slight deviation from this point will lead the system away, 
toward the stable fixed points. The stable fixed points here are for T > Tc (paramagnet: disordered spins) or T < Tc 
(ferromagnet: increasing ordering of spins). 

(ii) li < 0: hi and the components of hi decreases with A. Thus giving rise to irrelevant hi because larger A also 
implies larger correlation length. Here, the system moves toward a fixed point (attractive or stable fixed point). The 
fixed points at T = 0(< Tq) and T ~ oo(> Tc) stated above are stable fixed points. 

(iii) li = (vanishing eigenvalues): hi is called a marginal variable. In this case, one needs to expand Eq. ([2]) 
beyond the linear approximation (which was carried out earlier). Hence, Eq. ([7]) can be written as ^ 

A J^/ii(A) « Bh^iA), (9) 

where ztB (constant) leads to marginally stable {+B) and unstable fixed points {—B). For +B and A — !> oo, the 
solution behaves like Q], 

hAA) oc -^. (10) 

(iv) The eigenvectors that do not affect the physical properties or the components of hi are called redundant. In 
other words, redundant eigenvectors give only multiplicative constants. 

B. Coulomb screening: Shankar formalism 

We start by first giving a simple but very important example of rcnormalization, where the Thomas- Fermi screened 
Coulomb and the bare Coulomb potentials are actually renormalized potentials [3] . We follow the field theory formal- 
ism provided by Shankar [J], but the notations follow Ref. for easy comparison with the lET. In the field theory 
approach, the Coulomb potential, after introducing a smooth cutoff, exp[— |/c|/A] is given by [J] 



F(k,A) 



eo 



Ic^ -I- ^ K^ 



(11) 



where Kg is the Thomas-Fermi wavenumber and k is the wavenumber. It is straightforward to note here that 
VCk, 0) gives the Thomas-Fermi screened Coulomb potential, 

nk,0)= fj[j^,y (12) 



whereas, F(k, cxo) leads to the bare Coulomb potential, 



V{Koo) = —^. (13) 

Here, A is a cutoff parameter that has a dimension of k, and it is necessary to eliminate it from the final formula, 
either by taking A -^ oo or A — > as carried out above. The reason for doing so is that the final formula cannot be a 
function of the artificially introduced cutoff parameter, A. In the intermediate regime, k <C A -C (X^/k), Eq. PT|) is 
given by Q 

nk.A)=.-0j, (14) 

which describes the direct proportionality between T^(k, A) and A. This relation is very important to understand 
the formalism of lET with respect to RGT. 

Note 2: We will show how the lET can be used to capture the relationship between F(k, A) and A. A direct proof 
is given that formally associate the lET with the general formulation of RGT, as discussed earlier. The renormalized 
screened Coulomb potential, T4c plays the pivotal role with respect to strong e-e and e-ph interactions. 

C. Renormalization group equation: Gell-Mann Low formalism 

The GcU-Mann-Low equation is also known as the /3-function, or the differential RG transformation, is the principal 
result of the RG analysis Q. In fact, we have already introduced equations in the form of the Gell-Mann-Low 
equation many times earlier, which is given in Eqs. 1^, (jH]), ([3]), ([7]) and 0. Here, we will derive the equation in 
a straightforward manner. Assuming that we have a periodic lattice, an arbitrary coupling constant is expected to 
change or flow with respect to any changes to the cutoff parameter (A) that controls the distance between the lattice 
points (lattice spacing) or the number of lattice points considered (magnification) . For example, one can write the 
coupling constant (c) as Q 

c' = cA^-s^ (15) 

where g = 1 implies non-interacting case that leaves c unchanged. Furthermore, for g < 1 and g > 1 one have the 
repulsive (unstable) and attractive (stable) fixed points, respectively. This can be understood from the classification 
of eigenvalues discussed above: g < 1 gives larger c because 



.' ^ ^Al-g<l 



On the other hand, g > 1 gives smaller c due to 



A^-s<\ (16) 



oc cAi-s<o. (17) 



Usually, we need to set a cutoff value such that A = A,„in = or 1 or any finite values. Nevertheless, we will discuss 
this with respect to the lET in the following section. One then can iterate the RG equation, [Eq. ([T5|] to explore the 
sequence, c', c", ■ • •, and so on. Each iteration implies a change in the distance between the lattice points (spacing) 
or the number of lattice points as discussed previously. For example, Eq. (jlSp can be written as [5| 

c' = cA^-e, 

In- =lnAi-^ 
c 

c' = ce'"^'"', (18) 

Now, the flow of the coupling constant, c with respect to A is simply the differential equation given by 



dc 

(l-g)/dlnA, (19) 



c 

where this equation is known as the Gell-Mann-Low equation, or also called the beta (/3) function. All the solutions 
obtained from the lET Haniiltonian through the ionization energy based Fermi-Dirac statistics satisfy equations 
identical to Eq. (PO)) . 

II. CLAIM AND THE PROOF 

Prom the general formulation, and the renornialized Coulomb screening potential derivation given earlier, we know 
that there is a relationship between the effective Hamiltonian or F(k, A) and the length scale (A) [2|^, though one 
should remove it by means of the appropriate limits [for example, see Eqs. (fT2|) and (|T3| ]. In any case, one can actually 
expose a subtle issue here, which is: for a given system, the T^(k, A) can be A-dcpendcnt. which makes it possible to 
evaluate the evolution of the Coulomb potential between two different systems. For example, substitutional doping 
of element C into AB (ABi-ajC^r) implies that ^(k, A) should either be (i) dependent on A explicitly (provided that 
the other microscopic variables such as wavefunctions and Kg are fixed as constants) or (ii) the Coulomb screening is 
independent of A. If it is independent of A, then we need to know some details of the microscopic variables explicitly 
for each doping x. 

In lET, we have opted for the former strategy [(i) A-dependent] , which will be proven here. However, in lET 
formalism, we will not employ the field-theoretic cutoff-parameter approach since A is not only artificial but also 
microscopically not useful because it does not carry any information about the microscopic details of a particular 
system. For example, we will never know how to change A with increasing x, or for different doping elements. 

A. Claim 

From our earlier work, the many-electron atomic Hamiltonian is given by [6| 

' Z 



2to ^ ' 47reo ^ 



2|r,-r,| 



c{ri+rj) 



(21) 



where m is the electron's mass, e and Z are the electron's charge and atomic number, respectively. Moreover, 
I"! — r2| = \/rl + r2 — 2rir2 cos(02), while the labels, i and j identify the electrons and avoid counting them twice. 
The last term, which describes the screened e-e potential makes the Hamiltonian translationally not invariant for 
solids. However, the Hamiltonian can be made translationally invariant by assuming the dopants to be distributed 
homogeneously. For non-homogeneous distribution, we will make use of Eq. (j23p (Schrodinger equation with redefined 
eigenvalue) discussed below where we will leave Eq. (PT|) for atomic system, as it should be. The a is given by 0, 0] 

a = ;ie5^(-«'. (22) 

Here, n is the screening constant of proportionality, while A = (127reo/e^)aB, where eg Etnd Cb denote permittivity 
of free space and Bohr radius, respectively. In solids however, the many-body Hamiltonian can be written as [6| 

Hip = [Eo ± 0^. (23) 

In this case, we can prove for atomic system, Eq. (^5)) is equivalent to Eq. (PT|) @. Here, ^ is the energy needed for a 
particle to overcome the bound state and the potential that surrounds it. Eq denotes the total energy at temperature, 
T = 0. We claim that every solution obtained from the lET Hamiltonian [Eq. pT|) ] is a renormalized solution, which 
can be recast into RG differential equation, in the form of the Gell-Mann-Low formula. In order to prove this, we first 



need to show that the many-electron Hamiltonian given in Eq. ((2T|) satisfies the Zinn- Justin formaUsm. Secondly, one 
also needs to prove that the Shankar Coulomb screening formula given in Eq. p^ satisfies Eq. (P^ . Once we have 
done these, then it is straightforward to follow the proof of the claim. 

l>sc = ^e-'^-*''"'' = -j^e-^'. (24) 

B. Proof 

Let us recall Eqs. ([6]) and ([8]), which can be rewritten as [using Eq. ([23l) ] 

i 

h{A)y = Hip. (25) 

Here, AH\ is the slightly deviated Hamiltonian, from the fixed point Hamiltonian (H*) via linearization procedure 
given earlier. In addition, y and if are the eigenvectors, while A is known as the cutoff parameter that defines the 
RG iterations. We are trying to rewrite the general RG Hamiltonian of Zinn-Justin with the many-electron lET 
Hamiltonian [Eq. (HU]. In this lET Hamiltonian, the screened e-e interaction Coulomb potential is the one that has 
been rcnormalized and therefore, Eq. (|25p can be written as 

h{A)y = Vscip, 

h{A) - A'hil) ^ V,,{A) = A'Knsc, (26) 

where I is the eigenvalue (/ < 0: because Vqusc nrust decrease with A') that has its origin in the discreet linear 
operator, C* as described in Eq. (|3]). Here, Vunsc is the unscreened (bare) Coulomb potential that gives the strongest 
e-e interaction. Comparing Eq. ((26)) with Eq. (|24)) leads to the relation 

A?i-„ = A'^e-^-"*'^ (27) 

Now, taking the limit, I -^ gives A' — > 1, which in turn implies ^ — > cxd that leads Eq. ([26| to 

Fsc(A) ^ Knsc, (28) 

which is correct based on the lET. On the other hand, the limit, / -^ ^max gives A' — > A,„i,i (because / < 0) that 
implies ^ ^' and therefore we obtain 

V;c(A) ^ A'--V^^sc = Knsce-^^ (29) 

where, this result is also in agreement with the lET. However, for / -^ — oo gives A' — > that requires /i ^- oo and 
thus Eq. ([^ is given by 

V;c(A) ^ 0, (30) 

in which, this result satisfies both lET and the original result of the Thomas-Fermi screening. 

What we have done thus far is that we have proven that the cutoff parameter, A from the Zinn-Justin formalism can 
be directly related to the lET based screened Coulomb potential as given in Eq. ([27]) . In other words, V^c(A) = Kinsc is 
obtained for / = (A' = 1) that implies strong correlation between electrons. On the other hand, V|3c(A) = Kinsce"^'" 
is obtained for / = /,„ax (A' = Amjn) and subsequently, V^c(A) = for Z — > — cx) (A' — > 0). Therefore, the screened 
Coulomb potential, V^c is indeed a renormalized potential, without infinities. Infinities do not arise here due to the 
definition of Eq. (j27p . Note here that the limit, ^ — >■ oo implies finite value for Vunsc, which in turn gives strictly finite 



For example, KinscG"^'' can be related to the Thomas-Fermi screening potential by 
identifying Z^ax = ^tf and Ajf^^^ = Atf = e~^^ ■ Hence, from the lET formalism, we can find the evolution of the 
cutoff parameter, A for different doping, x and doping elements. 

Now one can proceed to obtain the Gell-Mann-Low equation. Using Eqs. ([^5]) and (P7)l . one can write 



Iny = lnA' = 



-fire 



4>? 



(31) 



dlnV 1, _l^. 



(32) 



d^ 



dhiV A/ire 2^^ 
where ln[T/sc/Vunsc] = \nV. Since InF = InA', we can rewrite Eq. (l33l) as 



(33) 



dC 



din A' 



d^ 



A' 



d^ 
dk^ 



Xfire 



iAC 



m), 



^'dXr-^(^)- 



(34) 



Equation (|34| has exactly the same form as Eq. ([20]) . 

The next step is to prove the lET version of the screened Coulomb potential, t^c can be transformed into the form 
of the Shankar Coulomb screening formula [Eq. ([T4|) ]. The Fourier transformed t^c [Eq- ([24]) ] is given by [3| 



Subsequently, one obtains 



v;c(k) 



or 



14c (k) 



eo 



'X + K^e-^i) 



(35) 



eoe 



-AC 



K? 



Fsc(k) 



eoAl 



eoK^ ■ 



(36) 



where, the limit k^ <C iiT^e '''^ has been invoked. Therefore, one obtains the relation. 



AiET = e^«, 
from Eq. (|M)) . Comparison between Eq. (|36p and Eq. P^ gives 



(37) 



Asha 



ikar 



kA 



lET. 



(38) 



Hence, Eq. ((38|) is the relation between the Shankar cutoff parameter and the lET-based cutoff parameter, where 
AiET is a function of the ionization energy or the energy-level difference (^) of the constituent atoms in strongly 
correlated matter. Importantly, when A — >■ 00, then we have Vj3c(k) — >■ 00, which is also in accordance with the limit, 
^ — > 00, as it should be. This divergence, 14c(k) -^ 00 is entirely due to the inequalities, k^ <^ K'^e~^^ (for Ajet) 
and k <C Ashankar ^ (-f'^s/k). Having found Eq. (P5|) . one can go on to prove that every solution of Eq. ([^51) is a 



renormalizcd solution. In fact, one just has to repeat the proof given below for every physical variable that has this 
exponential factor. However, we will only work out the proof for the V^c here, where this procedure can be easily 
carried out for other physical variables. 

Now one can proceed to obtain the Gcll-Mann-Low equation as follows: Eq. (|24p can also be written as 



AneQr 



ln^.-,.e-|^^.^, (40) 



' unsc 



where Eqs. ([55]) and ([5^ give us the procedure to recover the exponential function with the same mathematical 
structure as the Shankar smooth cutoff [J], which is given by exp[— |fc|/A]. Note however, Aiet is always dimensionless 
and refers to the distance between energy levels, unlike Aghankar that has wavenumber dimension and is based on the 
distance between the lattice points. Finally, one can derive the Gcll-Mann-Low differential equation from Eq. (PH)) as 
given below 

d\nV 1 ^ 1 ^f: ,,^, 

-^ . -A^re-.^^ (41) 



d\nV dV Xfire-i^^ 

Equation (|42p has exactly the same form as Eq. p4p . Now, we can repeat the procedure between Eq. p9p and 
Eq. P^ to obtain the Gell-Mann-Low differential equation for various physical variables mentioned previouslyB 

Equation p2|) also tells us that the flow of ^ is determined by the changes in the physical parameter, V. From the 
lET, we know how the changes in the elemental composition change £, a priori, and therefore, we can actually capture 
the consequential changes to the V. In other words, we do not need to know the changes occurring to V^ as a result 
of changes in the elemental composition a priori, in order to find ^. If we can find the systematic changes to V or 
any physical parameters mentioned earlier without knowing ^, then we do not need lET. 

Finally, from Eq. P^ . InV^ can only be independent of ^ if ^ — !• or ^ ^- oo. Therefore, one can obtain the 
corresponding /3-functions that satisfy these two limits (^ — > and ^ — > oo) as given below, 

2 

lim B(£t) — - — : satisfies Fermi gas, (43) 

lim /3(^J) — !• oo : superinsulator. (44) 

As a consequence, strong electronic correlation is captured by /3(^ > £,1) > 2/A/ir. 

III. DRESSED PHONON AND TOTAL DIELECTRIC CONSTANT 

The derivations covered in this section will be the backbone for the derivation of the acoustic and optical branches 
in the following section, and a further proof for the dressed phonon is given after that. Physically, £ in Eq. (|23p is 
the energy needed to excite a particular electron to a finite r, not necessarily r — > oo. In the early stages of the 
ionization energy theory, we used the atomic ionization energy (r — > oo) as the input parameter to compute carrier 
concentration [8J, and for this reason, ^ was labeled as the ionization energy. Furthermore, Eq. (j23p is technically 
easier to use, for example, since ^ is unique for each atom and using the total energy, Eo±£^ as one of the restrictive 
condition, we can derive the ionization energy based Fermi-Dirac statistics, which can be applied to solids. The 
respective distributions for electron and hole are given by Q 



where fcs is the Boltzmann constant. The dressed phonon frequency is given by [9[ 

where fip is the plasma frequency and K^ = 3nQe^/2eoEp. Here, Ep is the Fermi level at T = 0, fc and Kg are the 
wavenumber and Thomas- Fermi wavenumber, respectively Q- While, uq denotes the carrier density (free-electrons) 
at T = 0. The electronic dielectric constant based on the Thomas-Fermi approximation and ionization energy is given 

by0 

e(0,e,k) = H-^exp[A(i??,-e)]. (47) 

Using Eq. gll) and taking {K^Jk'^) exp[A(£'^ - ^)] > 1, we can rewrite Eq. (gll) as 






\k^~e%) 



(48) 



On the other hand, the constant of proportionality (e) that relates the Fourier transform of the total potential (0) 
in a metal to the Fourier transform of the external-charge potential ((/)oxt) can be written as [9| 

where, e is defined as the total dielectric constant, e(0,k). The ionic dressed dielectric constant, efo™'^''°'^ represents 
the screened ions due to screening electrons. We can rewrite Eq. (|151) . after taking (K^//s^) exp[A(i?j!. — ^)] ^ 1 and 
using Eq. (gT]) 

(t> = -0cxt == dressed — '?''cxt exp[A(£;^ - C)] , (50) 

In the presence of electronic medium with ionic medium as an external source, the potentials from the ionic ((/)ion) 
and the ionic induced electronic (0oxt) contributions can be written as [using Eq. (j47p ] 

eol0 = 0ion +0oxt, (51) 

Similarly, considering bare ions with electronic medium as an external source, we can write [9| 

ek)n™0 = </'cl + 0ext, (52) 

Noting that = 0oi + 0ext + (f'ion, and from Eqs. (|T7|) . P^ . (|5ip and ((5^ . we can arrive at 

,bare -i 

dressed -i i ion /i^Q^ 

'°" " 1 + (K2/fc2)eM^^-«) ' ^ ' 



If we take the ejjjj" = 1 — {flp/w'^), then it is straightforward to show that 

1 / 1 \ / w^ 



e 



1 + (X2/fc2)eMiJ^-0 ; \w^ - w{k) 



(54) 



Equations ((53|) and ((54)) suggest that in order to evaluate the changes of the strongly correlated electronic properties 
with different doping elements, Eq. (|Tf)) is sufficient. In this case, the electron-phonon interaction enters through the 
heavier effective mass [ll|. However, phonons cannot be treated as a mere variation to the effective mass if we are 
interested in thermal properties of correlated matter. 
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IV. DRESSED ACOUSTIC AND OPTICAL BRANCHES 



We consider a ID lattice with a basis with two different ions per primitive cell and in this system of diatomic linear 
chain, the harmonic potential energy can be written as [9| 



t/ha 



G 



^[ui(ria) - U2{na)]^ + —'^[u2{na) - Ui((n+ l)a)]^, 



(55) 



where ui(na) is the displacement of the first ion that oscillates about the site na, while for the second ion, the 
displacement is U2{na) that oscillates about na + d. Harmonic oscillation requires d < a/2. Using Eq. (|50)) we can 
write the interaction potential constants, Q and G as 



Q = 



G 



d^(j){na) 



exp[A(^ - E%)] 



d'^(p{na + d) 



exp[A(C - E%)] 



(56) 

(57) 



Here, ^ is the average ionization energy in the vicinity of ionic mass 1 (Mi) and 2 {M2). (f>{na) and (j){na + d) 
represent the interaction energies between the respective ions through the screening electrons. The exponential term 
originates from our dressed electron described earlier. The equation of motion based on Eq. (|55|) can be readily solved 
to obtain 



W4 



1 



2M1M2 
X (Afi + 



(Q + G)(Afi+Af2)± [{Q + G) 



M2Y 
2QGcos{ka) 



4MiM2[{Q + Gf - Q^ - G^ 
1/2-1 



(58) 



It is clear here that the ionization-energy dress that appears in Eq. ([55)1 is identical with Eq. pS)) . In the limit 
^ —7- 00, then the above system is infinitely rigid. In view of Eq. (^5)) and (|58p . we can readily derive (following 
the procedure described in Ref. [9|) the ID harmonic oscillator Hamiltonian (in the second quantized form) and its 
eigenvalues, respectively as 



H 



k 



?iw(k)e2 



hH^-K) 



aj^flk 



(59) 



i? = ^;ia;(k)e^^(«-^-) 



nk+2 



(60) 



Note here that ?7k is the dressed phonon distribution function. The a]^ and au denote the usual phonon creation 
and annihilation operators. 



V. DRESSED ELECTRON-PHONON INTERACTION 



In the previous sections, what we have done is that we have stripped the ionization-energy dress (exponential term) 
from the electrons and put it on the phonons (dressed phonon) . This implies that we can readily tack the exponential 
term to the phonon frequency and derive the Debye model and heat capacity. But before we move on, we need 
to prove that such discriminative dressing (stripping the electrons in order to dress the phonons) is mathematically 
and theoretically valid. We can prove this by showing that when we tack the exponential term (ionization-energy 
dress) onto the phonon frequency, we will find that the form appears in the derived dressed e-ph interaction potential 
is exactly the same as the form, appeared in the dressed electron-electron (e-e) screened Coulomb potential (given 
below) 
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yec(k) = -^ 



We know from the second-order time- independent perturbation theory [12|, 



(61) 



E 



(2) 









£;(0) 



£^, 



(0) 



(62) 



where i?"' is the second order correction to E^'^K The tfm,n denotes wavcfunction for states m and n, respectively, 
while -ffpcr is the perturbation. If we take i?por as the e-ph interaction potential, iJcp for solids in k-space, and using 
Eq. gH) then we can rewrite Eq. ([62]) 



S(2 



E 

k,k' 



|(</7k|i?op|Vk-) 



£:(k)-S(k*)-rK^(k-k*)e 



hK(.-E%) 



y^ ?^k(l -J^k*) 



k,k* 



|(<^k|-f^cp|<^k- 



S(k) - S(k*) - ritj(k - k*)e5^(«-^?-) 



(63) 



The detailed reasons why Eq. (|62p represents the e-p/i interaction and the missing first-order correction (its matrix 
elements equal zero due to orthogonality) are given in Ref. [lO|. The dressed phonon distribution function can be 
written as 



1 



nk 



o/3?icj(k)e2^ 



(64) 



The denominator in Eg. (j63p have made use of the crystal momentum conservation, k — k* = q, and /3 — l/kgT. 
Using Landau approach [9|, we can derive the e-ph interaction potential, Vk_k* from Eq. (j63p 



Vl 



k,k* 



_d^E^ 

9nk9nk* 



l.9k,k- 



2ncj(k-k*)e 



^A(C-BO) 



[huj{k - k*)e^^«-^^)]2 - [^(k) - E{k*)] 



The e-ph coupling constant, l^k.k* | = K'Pk |^op| Vk*)| is given by 



l5k,k- 



Veo{\k~k*\^ + K^)2 



hu{k^k*). 



(65) 



(66) 



In the previous work [11| . we have claimed that the ionization energy dressed electron-electron screened Coulomb 
potential has the electron-phonon interaction effect taken into account implicitly that gives rise to heavier effective 
mass. Consequently, we need to impose this condition of heavier effective mass into Eq. (|65|) . Here, we will make use 
of the effective mass {m*) theorem [9|, in which it is approximately given by 



d^Enik) 



oc 



o 



(67) 



dk,dkj 

This means that for heavier effective mass, we need to satisfy, [huj{k — k*)e2 -'^(i~^F)]2 ^ [£'(k) — E{k*)]'^. In other 
words, we need the term, [E{k) — E{k*)]'^ to be relatively small in order to comply with Eq. (|57)) . which gives rise 
to the heavier effective mass effect. As such, Eq. (|65|) can now be rewritten as (also, after dressing Eq. (|66)) with the 
ionization energy) 



Kp(k,k*).- 



Ik-k* 



K^ exp[X{E^p ~ 0] 



(68) 



As anticipated, Eq. (|68p is in exact form with Eq. (|6T|). This completes the proof that the ionization energy is 
proportional to the electron-phonon coupling constant and its relationship with heavier effective mass, where we have 
only explained this proportionality qualitatively in Ref. |ll| . The next step is to apply the outcome of this proof to 
derive the heat capacity and the well known Debye model (Debye frequency and temperature). 
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VI. DRESSED HEAT CAPACITY AND DEBYE TEMPERATURE 



Now, we will make use of all the results derived earlier to derive the dressed heat capacity formula, dressed Debye 
frequency and the dressed Debye temperature. The energy density of a harmonic crystal is given by (using Eq. (|60p 
and introducing the static potential energy constant) 



U =[/.tat + ^E^^«We^'^'"''^^ + 






k.s 6' 



(69) 



and the non-constant heat capacity (C„) is given by 



k 



d hujsCk^e^ 



5^(?--Ef) 



V ^ dT g/3Rcj,(k)e^- 



(70) 



We also have introduced the branch index, s — l,...,3p, where p is the number of ions in the basis signifying that 
for each k there are 3p normal modes 0. In the limit of long wavelength (small fc), w(k) = Cs(k)fce2'^^^^^i') [using 
Eq. (|48p] where Cs is the long wavelength phase velocity that depends on the vibrational polarization, s and the 
direction of k [^. Again, using the same procedure given in Ref. Q, we solve Eq. ([70)1 and obtain 



L^v — 



27r2fcf 



kuT 



l3 



he 



.-Im-K) 



(71) 



where fcs is the Boltzmann constant and note here that Eq. (j7ip correctly gives the inverse proportionality between 
Cy and the ionization energy, ^. Later, this proportionality will be associated to the stiffness of the crystal, as we can 
anticipate from this relationship itself. 

Next, let us work within Dcbye's approach and derive the Debye frequency, temperature and the heat capacity, all 
dressed with ionization energy. In this case, we may simply write. 



c 
where c is the constant sound velocity. The phonon density of states, N(uj) is given by 



(72) 



N{io) = — = — 
doj du! 



i:^ -iA(c-B°) 



Gtt^c-^ 



27r2c3' 



(73) 



where N is the number of phonons. By defining lu = lu-£, and rtph = N/V as phonon density, we can rewrite Eq. (j73p 
to obtain the Debye frequency 



WD ==c[67r^nph]^/'^e2 



l/3p|A(?-B«.) 



(74) 



Also noting that N/V = [k^e 2^(? ^f^^/Gtt^ and ksQu = f^D = hckj^, then we can arrive at the Debye 
temperature 



Od 



he 



[eTT^riph]^ 



ll/3piA(5-_E«,) 



(75) 



Before we continue, it is important to reduce the level of confusion that may arise in the above formulations and 
proofs with the exponential term appearing here and there. Hence, we will give a step-by-step derivation for the 
Debye heat capacity, which will prove the origin of the exponential term in Eq . (|75p , as well as points out the correct 
way of handling this exponential term. Again, we start from Eq. (|70p using [9| dk = k^dkdil to get (il is the phonon 
polarization volume in k-space) 
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Uy 






J^(i-E") 



Now we make the substitution, x = fiKcke'^^^^^^p^ into Eq. ([75)1 and obtain 



(76) 



a, = 9/s 



B?^ph 



T 



eDe^^(«-^^) 



3 rajn 4 a; 

' ^ X e 



(e- - 1)^ 



-da;. 



(77) 



We can see that Eq. (|77p is in agreement with Eq. ([TT]) and ([75)) . 



VII. CONCLUSIONS 

In conclusion, we have proved that the ionization energy theory is indeed a renormahzed theory based on the 
approach developed by Zinn- Justin, Shankar and Gcll-Mann— Low formalism. Subsequently, wc have applied these 
proofs to show that the ionization-energy dressed clcctron-phonon interaction potential is exactly in the same form 
as the ionization-energy dressed electron-electron Coulomb potential. Consequently, we can obtain the heat capacity 
and the Debye frequency as a function of ionization energy. These latter proofs can be used to explain the evolution 
of the electron-phonon interaction with respect to elemental composition in strongly correlated matter. Apart from 
that, this theory may define the existence of the (quantum) phase transition between free-electron metals and strongly 
correlated matter, which will be addressed in our future work. 
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VIII. APPENDIX 
A. Reversible transformation between Eq. (|46p and (|48[l 



One should not assume that Eq. (|46)) and Eq. (|48)) are not identical or theoretically different. Equations (|46l) and 
Eq. (|48p arc both mathematically and theoretically exact, which can be transformed and reverse-transformed exactly. 
We have to prove this transformation because obtaining Eq. (P5|) from Eq. pS)) is our starting axiom for the later 
proofs in this appendix. Hence, by taking the appropriate limits. 



hm e(0, e, k) = 1 + § cxp[A(i??, - i)\ . 
i>E% k^ 



lim e(0,e,k) = l + §=e(0,k). 



and also after taking, {K^/k'^) exp[A(i?^ - ^)] > 1 and K^/P > 1, we obtain 



lim a;(^,k) = ^ exp 



^KC-e'f) 



lim w(^,k) = M^ = ,^(k) = "^ 
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For conventional free-electron metals, the linij_^^o applies and any perturbation due to temperature and/or po- 
tential (due to doping, defects, etc.) will vary the Ep and uq accordingly. For example, for each type of free-electron 
metal, and for each type of perturbation, there arc uq and Ep, which arc both unique. In our approach, any variations 
to Ep and uq are captured through ^. Meaning, for a given system, Ep and np are constants, while ^ will respond to 
any perturbations, be it large or small. For example, n(r, T, doping) = no exp[A(£'^ — £,)] that will tell us the changes 
to n(r, T, doping). In other words, n(r, T, doping) = no (of that perturbed system), or simply, lini|^£;o n(r, T, doping) 
= no (of that perturbed system), which will take us back to square one, the free-electron metals Q- As a consequence, 
we can see that a;(^,k) — a;(k) or Eq. (|46| and Eq. ()48|) are exactly the same due to our reversible transformation. 
Simply put, the limt^^o implies Ep and no have been taken as constants for a given system and we are letting ^ to 
respond accordingly to any external and/or internal perturbations. 

B. Derivation of Eq. (f54)l 

We start from the definition as given in Eq. (|49| where, 

e</' = </'cxt. (78) 



Next, we add Eq. dST]) and Eq. (|52]) to obtain, 

0(eel + ^oT) = 20ext + </'ion + (j^el- (79) 

We then substract Eq. ^ from Eq. dTQ]) to get 



(/.(eel + e^r - e) = '/'oxt + <^io„ + <^ci. (80) 

Since, 4> = 0oxt + (f'ion + 4>c\-, Eq- (|80| needs to satisfy 

e = eei + ej;r-l- (81) 

Now, if the electrons screen the dressed ions, then Eq. ([75)) can also be written as 



such that [from Eq. ^ and Eq. ([ST]) ] 



efoT'^'eeO^ = </'ext, (82) 



dressed i_ bare -i (il'X\ 



Therefore, 



,barc 1 

1 + J2S . (84) 

Eel 

Substituting Eq. (gT]) into Eq. dHU) will lead us directly to Eq. ^^. Subsequently, we substitute ej5,n" == 1 ~ (^p/'^^) 
and Eq. (|i5)) into Eq. (jM)) so as to arrive at 

^bare i 



Col 



Cel 
= 1 



(1^2/^2 



Col 

2, 



_ ^ (»^/ecl) 



1-^- (85) 
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Prom Eq. (|83|) . we can rewrite the total dielectric, e as 



^ ^drcssed^ 



(86) 



Substituting Eqs. ([53]) and (gT]) into Eq. dHU) will lead us to Eq. dSH) as shown below, 

1 1 1 



1 



1 + (x2/fc2)e^(^5^-«) )W- «i(k) 

C. Derivation of Eq. ([5g)l 

Using Eq. (j55p , the equations of motion can be written as 

5f/harm 



(87) 



MiUi{na) = — 



dui{na) 

Q .r 

2 ^ ^ 
= — Q[wi(na) — U2{na)] — G'[u2(na) — wi((n + l)a)]. 



— ■ 2[ui{na) - U2(na)] - — • 2[u2{na) - Ui((n + l)a)] 



Ai2^t2(?^a) 



9^t2(f^a) 

O G 

= -— • 2[ui{na) - U2{na)\{-1) - — • 2[u2{na) ~ ui{{n + l)a)] 

= —Q[u2{na) — ui{na)] — G[u2{na) — ui{{n + l)a)]. (89) 

The dressed harmonic potential energy (0'^'^°*'®°'^(a;)) and the interaction potential constants (Q and G) are respec- 
tively, given by 

^drcsscd(^) = <^(x) exp[A(C - E°p)], (90) 



Qd.csscd ^ ^'^"[(^«)'^] = Qexp[A(e - EU, (91) 

G'i--<='i = a'0'^"'^'^°'^[(»Q + rf),x] ^ g^^pj_^(^ _ ^0 )]_ (92) 

Here, the two ions with their respective masses. Mi and M2 have been considered, in which the ionic oscillations 
are harmonic (d < a/2). The solutions we need are in the form of 

ui{na) = fo^e*^'^""-"*), U2{na) = 62e'(''"°""*\ (93) 

where k is the wavevector, while 61 and &2 are constants. By substituting Eq. ([M|) into Eq. ([55)1 and (|89p . one can 
obtain two coupled equations: 

[MiLj^ - (g + G)e^(«-^^)]&i + [g + Ge''^'^]&2e^(«-^^) = 0, (94) 



[Q + Ge^'^^lbae^^?-^^) + [A/sw' - (Q + G)e^(«-^^)]6i = 0. 
These homegeneous equations can be solved by setting the determinant equals zero, i.e., 



[Nhio^ - (Q + G)e^(«-^-)][Af2C^' - (0 + G)e^(«-^5^)] = [(Q + Ge''='')62e^(«-^^)]^ 
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(95) 



(96) 



Therefore, 



[MiM2]w4 + [-(Q + G)(Mi + M2)e^(«~^^ V 

+ [(Q + G)2 - Q2 _ ^2 _ 2QGcos(fca)]e2^(«-^^) = 0. 



(97) 



By identifying A = AhA'h, B = -(Q + G)(A/i+M2)e^(«-^S-) and C = {Q + Gf -Q^ -C^ -2QGcos{ka)y^'^^~'^'^\ 
we can obtain Eq. ([58| from Eq. (|98|) given below: 



B ± \/b2 - 4AC 



2A 



(98) 



D. Derivation of Eq. ([68|l 

From Eq. (p^ . -E,i and Em are the zeroth order initial and excited states, respectively. As such, their difference 
due to phonon absorption is given by 



ii;(°) - i?(°) = E{\.*) + hLo{), k*) - i?(k) 



(99) 



Prior to phonon absorption, one can identify the excited states as unoccupied, while the initial states as occupied, 
thus, we need to include a probability factor in the form of X^k k* ^k(l — '^k*), which satisfies the Bosc-Einstein 
statistics. As a consequence, we can arrive at Eq. (p5)) . The next step is to obtain Eq. (p5|) : 



Fw w = 



92 






■ ^ nk(l - 71k* 

k,k* 

92 



Kvk |i?op|<pk*) 



^TikaTik- ,^. '""'^ ■"" £^(k) - £;(k*) - huj{\i - k*) 

7ik(i-"-k*) , 'ik*(i-"k) 



9nk' 
2 



9nk9nk* 

|2 ^ 

-|5k,k- 

Iffk.k*!^ 



£:(k)-£;(k*)-fiw(k-k*) £;(k*) - £:(k) - fii^(k* - k) 

(l-nk-) TT-k* 



£:(k)-£;(k*)-fe^(k-k*) s(k*)-s(k)-ntj(k*-k) 
1 1 



£:(k)-s(k*)-rK^(k-k*) s(k*)-s(k)-ncj(k*-k) 

2fi^(k-k*)e^^(«-^-) 



. [;iCj(k - k*)e3M?-BS.)]2 _ [£;(k) _ ^(k*)]2. 

where |(?k,k'|^ = Iffk'.kP- After dressing Eq. (|66p with ionization energy, one obtains 



1 



'^•''''* ' V eo [|k - k* |2 + if 2gA(5-£0 )] 2 

Substituting Eq. ([TOT|) into Eq. (fTOO| leads to 



ific^(k-k*)e5A(«-B^). 



(100) 



(101) 
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K=p(k,k*) 



2V^eo [|k-k*|2 + if2eA(C-i=;°)] 



2fii^(k-k*)e5M?-^^) 



[;iw(k - k*)e^^(«--^?-)]2 - [£;(k) - ^(k*)]2 



(102) 



Now, invoking the effective mass theorem [Eq. ([GT])]. Eq. (|102p can be written as given below to obtain Eq. ([68 



Kp(k,k*) 



2Feo [|k-k*|2 + X2e^(«-s^)] 



2;iw(k-k*)e5^(«-^^) 



[;ia;(k - k*)e^^(«-^^)]2 



1 



|k-k*|2 + /^|exp[A(S0 -0] 



(103) 



E. Derivation of Eq. ([71]) 



Equation (|5(I)) can also be written as 



U =C/,,,, + l^lfi^,(k)e^^«-^-) 



y ^^ 2 



l^nc.,(k)e^^«-^^)nk, 



(104) 



After substituting Eq. (|M)) into Eq. (|104p . one can arrive at Eq. (|M)) . Whereas, the C^ can be derived from Eq. (|104p 



_aC/_j_^ (9 fiij,(k)e5M«-EF) 



k.s ^ J^ 

^ d f dk ;iw,(k)e3^(«-^S-) 

2^ AT 



9T J (27r)3 g/3Rc.,(k)e^^'«-^F' _ i' 

Next, we substitute these equations given below into Eq. (|105p : 



(105) 



;(k) = Cs(k)A:e^^(«-^^\ dk = k^dk sin 6* dO dcj) = k'^dkdn, 



(106) 



dk 



(107) 



Hence, we get 



C'y 



d f Pdk dn ;ic^fce^^(«-^^) 



kMi-E") 



^ dT J (27r)3 g/3rxc,fee^"- -f' _ ^ 

dn r k^dk fic,/ce5^(«-^S-) 
47r 






27r2 g/SRc.fcei^K-^?-) 



1 



—y-f 

27r2 ^ dT J 



dn 

47r 



x^dx 



_3__a_ 

2^aT 



^E 



[/3ft, 
df} 1 
47r c.< 



c,e^^(«-^^)f 



x//? 



/34fi3eiA(C-£0) 



oo ^3 



X dx 



(108) 



Using 
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I 1^ r dn 1 



(109) 



/•°° x^dx _^ r 



3"-"^da; := 6 



X e 



n=l 



TT 

l5' 



we obtain Eq. ((7T|) . 



O^j 



TT^fcs ar* 



10 ar 



/ic 



-|A(C-B?.) 



27r2fcf 



Alternatively, instead of using Eq. (|107p . we use 



kr,T 



he 



-|A(«-_E«) 



e 2 



dx 



X = phcgk, — - = /3/iCs 



to get [also after using Eq. (|109p ] 



_ d .i-^ f dn f Pdk hc,ke^^^^-^'^'> 



ar 



47r 



27r2 ^/3;ic,fccxp[iA(5-_E«)J _ j^ 



1 d ^ f dn 

3 e^^(«-^^) 9(l//34) 





2 


dx 
[l3hcs\ 



xe^ 



kHi-Ei-) 



27r2 ft3c 



dT 






xp [iA(C- 



■£«: 



Substituting 



(110) 



(111) 



(112) 



(113) 



^oo 3 7 oo „oo r -1 -1 4 oo 



1 
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into Eq. (J113p . we can see that Eq. (|113l) becomes Eq. piip . which is Eq. (fTTI) . 



(114) 



F. Derivation of Eq. ([75]) 

From the left-hand side (LHS) equation of Eq. (|106p , 



k = —e^^^ii^Ep)^ 



On the other hand, the number of phonon modes, N is given by 



(115) 



N ■■ 



(47r/3)fc3 _ Vk^ 



(27r/L)3 67r2 ' 



(116) 



19 



where L and V arc the length and volume, respectively. Next, by substituting Eq. (|115p into Eq. (J116p . 



Ar.Z^e-i^(«--^), (117) 

and differentiating Eq. (|117p with respect to ui will lead us directly to Eq. ([75]) . Substituting Eq. (|115p into Eq. (|117p 
gives 



Subsequently, we define, k = kjj and w = a;^), thus 






(118) 



ks^D = huj£i ~ hcki), Qd = - — kjj. 

ks 



Inserting Eq. (fTTS)) into Eq. (fTTO)) gives 



ks 



6^2^g|A(C-B°) 



V 



(119) 



(120) 



which is equivalent to Eq. ((75)) by noting that riph ~ N/V. 

G. Derivation of Eq. ([77)) 

From Eq. dTO]), noting limy_>oo(l/V^) J2k = I '^k/(27r)3, and using the RHS equation of Eq. ([1M| one gets 

1 



where 



a =p£r%3^fc^iA(,-.=^)5 



27r2 



9r 



g/9Rcfee5^<«"-^F: 



1 



(121) 



^E 



dJ7 

47r 



- isotropic 



cs(k) 



(122) 



polarization 



_d_ 
df 



1 



^fihckeh^^i-^%) 



ar' 



^iShcke^ 



h>-{i-E%) -.-I 



1]" 



-![' 



,,3ncfce7 '"-''- 



^(«--Ef) -1-2 



^'^^ /3?tcfce^^"^~-^f ' 



^icfc e'^'^' 



,cfeei^'«-^F' 



fcBT2 



Therefore, we can derive Ct, as given in Eq. (|76p . Consequently, after the change of variables as given in Eqs. 
and dMl), 



ksT^ re/3ftcfcei^<«-^F' _il' 



(123) 
(fTMll 



C„ — 



3(;ic)2 /■''" ,, „ ,^f_B.o^ e- 



fc4dfce^(«-^-)- 






2TT^kBT^ Jo [g/3?icfce7M5-i=^) _ 

3{hc)^ r° (a;fcBr)4e^(?-^F) dxfcsT e 



Phcke^'- 



2n^kBT^ Jo (nc)4e2^(«-^^) fce^^^?"^^) [g/3ncfce5^«-^^' „ i]2 



3(;ic)2 


'keT' 
he 


2Tr^ kgT^ 



n5 



(l-2-i)A(e-B5!,) / ° ^ g^rfx 



e^ - 1 



(124) 
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Prom Eq. (J119p . one gets he = {ksOD/kD)- Substituting this and Eq. (|116p into Eq. (|124p . and after some algebraic 
rearrangements, one can arrive at 



a 



3{hc) 



27r2/cBT2 
27r2(nc)3 



he 



,(l-2-i)A(C-B«,) / ^ g 



(e- - 1)^ 



B± g-fA(C-B^) 



(e^ - if 



^(6^'nph) 



T 



enei^(«-^^) 



^ '■"° x^e^rfx 



e^ - 1 



,2' 



(125) 



which is nothing but Eq. ^7^. Again, using Eq. piOp and xd = Phcko = [kB^o/kBT) = Qd/T -^ cx) for T < Qd, 
one can solve Eq. (|125p : 



(e^-iy 



E 



Jo ±^.n* 



247r^ 
90 ' 



(126) 



therefore, 



C^ = — r— nphfcs 



T 



l3 



ez3pi^(«-^s^) 



(127) 



If we use Eq. (J112I) instead, then we need to invoke Eq. pi4p and the result [Eq. (|127p ] remains the same: From 
Eq. (USD 



a 



_3{hcf_ 
2TT^kBT^ Jo 



k^dke^^-^'^^^ 



3(M' xu-E%) 



2t:^ ksT^ 



gl^hcke'! 

4 



^a-E") 



Phc 



dx 



Phc 



-IF 

g^A(«-EO.) 



^^eiAK-Bj,) _ _^j2 



3fc^ T'^ " ''^° ^4^a;e5^(« -^F' 
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we can rewrite Eq. p29p as 
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which is Eq. (fT27)) . 
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